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Abstract 

Although the spectra of heavy quarkonium systems have been successfully 
explained by certain QCD motivated potential models, their strong decays 
are still an open problem. We perform a microscopic calculation of vector 
charmonium strong decays into open-charm mesons where the qq pairs are 
created from the same interquark interactions acting in the quark model that 
has been used to describe its spectrum, and also its leptonic and radiative 
decays. We compare the numerical results with those predicted by the ^Pq 
decay model and with the available experimental data, and discuss the pos- 
sible influence on the strong widths of the different terms of the potential. A 
comparison with other predictions from similar microscopic decay models is 
also included. 
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1. Introduction 

Since its discovery in 1974 [ll, [ij, the charmonium system has become the 
prototypical 'hydrogen atom' of meson spectroscopy js], 0]- The spectrum 
of relatively narrow states below the open-charm threshold at 3.73 GeV can 
be identified with the IS*, IP and 25* cc levels predicted by potential mod- 
els which incorporate a color Coulomb term at short distances and a color 
confining term at large distances. Difficulties arises above the open-charm 
threshold due to the coupling with the meson-meson continuum. 

In recent times the construction of 5-factories has triggered the discovery 
of many new particles, denoted as XYZ mesons, whose nature supposes a 
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challenge for the theorists and can contribute to a better knowledge of the 
cc phenomenology. 

One open topic of great interest on charmonium states is their strong 
decays which constitute a rather poorly understood area of hadronic physics. 
A great part of our knowledge of strong interaction comes from decays and 
for that reason it is important to pursue the most complete description of 
them. 

Attempts of modeling strong decays date from Micu's suggestion \E] that 
hadron decays proceed through qq pair production with vacuum quantum 
numbers, J^*" = 0^+. Since this corresponds to a ^Pq qq state, it is now 
generally referred to as the ^Pq decay model. This suggestion was developed 
and applied extensively by Le Yaouanc et al. 0, 0] in the 1970s. Studies of 
hadron decays using the ^Pq model have been concerned almost exclusively 
with numerical predictions, and have not led to any fundamental modifica- 
tions to the original model. Recent studies have considered changes in the 
spatial dependence of the pair production amplitude as a function of quark 
coordinates [8-ll|. There have been some studies of the decay mechanism 
which consider an alternative phenomenological model in which the qq pair 
is produced with ^Si quantum numbers 12|. However, this possibility seems 
to disagree with experiment [9|. 

An alternative procedure is the study of strong decays through micro- 
scopic decay models. The difference between this approach and those de- 
scribed above lies on the description of the qq pair creation vertex. In the 
microscopic decay models the qq pair is created from interquark interactions 
acting in the quark model. The differences between calculations of this kind 
lies in the choice of the pieces of the potential which enter in the vertex 
calculation. 

So Eichten et al. 13|, assumes that the qq production is due to the 
time-like part of the vector Lorentz confining interaction, while Ackleh et 



al. [15(1 and Bao et al. [16|| assumes that the qq pair is produced by one- 



gluon exchange and scalar confining interactions. Then, the description of the 
strong decays is intimately related with the problem of the Dirac structure 
of the confinement. 

From the point of view of the hadronic spectroscopy, confinement has 
to be dominantly scalar in order to reproduce the hyperfine splittings ob- 
served in heavy quarkonium although to explain the light quark 
phenomen olog y a small mixture (of the order of 20%) of vector confinement 



is needed [20|. However, build Hamiltonian-based models of QCD |21 



seems 
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to require vector confinement. 

In the present work, we generalize the microscopic decay models of ref- 



erences mentioned above using the quark interaction of Ref. [20[ which in- 
cludes one gluon exchange plus a mixture of scalar and vector confinement. 
This model successfully describes hadron phenomenology and hadronic re- 
actions and has recently been applied to mesons containing heavy quarks in 



Refs. [22|, |23[. We study the open-charm strong decays of the 1 cc reso- 
nances looking for the possible influence of the mixture of scalar and vector 
Lorentz structure. 

The content of the present paper is organized as follows. In section [2] we 
briefly explain the main features of the constituent quark model. Section [3] 
is devoted to the microscopic description of the strong decay mechanism. 
Section H] is dedicated to show our results and how its compared with those 
of the ^Pq model. A comparison of the numerical values coming from our 
microscopic decay model with those of other similar microscopic decay models 
is also included. We finish in Section |5] with some remarks and conclusions. 

2. Constituent quark model 

One consequence of the spontaneous chiral symmetry breaking is that 
the nearly massless "current" light quarks acquire a dynamical, momentum- 
dependent mass M{p) with M(0) ~ 300 MeV for the u and d quarks, namely, 
the constituent mass. To preserve chiral invariance of the QCD Lagrangian 
new interaction terms, given by Goldstone-boson exchanges, should appear 
between constituent quarks. This together with the perturbative one-gluon 
exchange (OGE) and the nonperturbative confining interactions are the main 



pieces of our potential model [20 



The wide energy range covered by a consistent description of light, strange 
and heavy mesons requires an effective scale-dependent strong coupling con- 



stant. We use the frozen coupling constant of Ref. |20 



In(^) 

where /i is the reduced mass of the qq pair and ao, fio and Aq are parameters 
of the model determined by a global fit to all meson spectrum. 

In the heavy quark sector chiral symmetry is explicitly broken and Goldstone- 
boson exchanges do not appear. Only the OGE and confinement potentials 
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arc present and contain central, tensor and spin-orbit contributions. For the 
OGE potential, they are given by 




(2) 



where S± = ^(c?i it (Xj). Besides, ro(/i) = and rg{fj,) = are 

regulators which depend on Hij, the reduced mass of the interacting quarks. 
The contact term of the central potential of one-gluon exchange has been 
regularized in a suitable way as 



The breaking of the color electric string between two static color sources 
is a phenomenon predicted by QCD and it is the basis of the meson decays 
and hadronization processes. Although there is no analytical proof, it is a 
general belief that confinement emerges from the force between the gluon 
color charges. When two quarks are separated, due to the non-Abelian char- 
acter of the theory, the gluon fields self-interact forming color strings which 
bring the quarks together. 

In a pure gluon gauge theory the potential energy of the qq pair grows lin- 
early with the quark-antiquark distance. However, in full QCD the presence 
of sea quarks soften the linear potential, due to the screening of the color 
charges, and eventually leads to the breaking of the string. We incorporate 
it in our confinement potential as 



rij/ro 



S{fij) 



(3) 




[-a,(l-e-^=^-) + A] (A^Ap, 




(4) 



+4mim,(l - as)){S+ ■ L) + (mj - ml){l - 2a,)(5_ • L) , 
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where controls the ratio between the scalar and vector Lorentz structure 

VcoNir-,) = asVS'S^^^in,) + (1 - a,)V^^oT ir^,) . (5) 

At short distances this potential presents a linear behavior with an effec- 
tive confinement strength a = — Oc /ic (A^ • Ap while it becomes constant at 
large distances. This type of potential shows a threshold defined by 

^^/.r = {-ae + A}(A^Ap. (6) 

No qq bound states can be found for energies higher than this threshold 
and the system suffers a transition from a color string configuration between 
two static color sources into a pair of static mesons due to the breaking of 
the color fiux-tube and the most favored subsequent decay into hadrons. 

To find the quark-antiquark bound states with these interactions, we 



solve the Schrodinger equation using the Gaussian Expansion Method [24 
It allows us to evaluate easily the strong decay amplitudes. 

In this method the radial wave function, solution of the Schrodinger equa- 
tion, is expanded in terms of basis functions as 

Rair) = <€i{r), (7) 

n=l 

where a refers to the channel quantum numbers. The coefficients and the 
eigenenergy E are determined from the Rayleigh-Ritz variational principle 



E 

n=l 



(8) 



where T",„, N^,^ and V^"' are the matrix elements of the kinetic energy, the 
normalization and the potential, respectively. T°,^ and N^>n diagonal 
whereas the mixin g be tween different channels is given by V^'^ . 



Following Ref. [2J] we employ Gaussian trial functions whose ranges are 
in geometric progression. This is useful in optimizing the ranges with a small 
number of free parameters. Moreover, this distribution of range parameters 
is dense at small ranges which is well suited for making the wave function 
correlated with short range potentials. The fast damping of the gaussian tail 
is not a real problem since we can choose the maximal range much larger 
than the hadronic size. 
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Quark mass 


rric (MeVj 


1 'IPO 


(iOTl Tl Tl PTTl P'^"^ 


n (MpV) 








0.576 




A (MeV) 


184.432 




as 


0.81 


OGE 


«o 


2.118 




Ao (fm-^) 


0.113 




yuo (MeV) 


36.976 




fo (fm) 


0.181 




rg (fm) 


0.259 



Table 1: Quark model parameters. 



J^^' n MThc. Mexp. 



1 


3096 


3096.916 ±0.011 


[25] 


2 


3703 


3686.108+|]:[5}^ 


[25] 


3 


3796 


3778.1 ± 1.2 




4 


4097 


4039 ± 1 




5 


4153 


4153 ±3 


[25] 


6 


4389 


4361 ± 9 ± 9 


[26] 


7 


4426 


4421 ± 4 


[25J 


8 


4614 


4634+?^^ 


m 


9 


4641 


4664 ± 11 ±5 


[26] 



Table 2: Masses, in MeV, of 1 charmonium states. We compare with the well established 
states in Ref. [25[ and assign possible XYZ mesons. 



Table [T] shows the model parameters fitted over all spectrum of mesons 
and needed for the heavy quark sector. Note that, once the parameters are 
fitted, our confinement interaction is dominantly scalar. 

Finally, Table |2] shows the masses predicted by our model for the vector 
charmonium states, the comparison with the experimental data and some 
possible assignments of XY Z mesons. Further details on the spectrum and 



other properties of vector charmonium states can be found in Ref. [22 
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3. Strong decays 



The microscopic decay models are an attempt to describe strong interac- 
tions in terms of quark and gluon degrees of freedom. As mentioned above, 
after the pioneering work of Eichten et al. 13|, which assume that the 



strong decays are driven by the time-hke component of the confining inter- 
action, only few works have addressed this topic in a partial way without 
discussing the relationship between this decays and the nature of confine- 
ment. 

We shall assume that the responsible for the strong decays is the full 
quark-quark interaction of our model which includes one-gluon exchange, 
scalar and vector confining interactions, allowing in this way the study of the 
influence of the different pieces on the final results. 

The associated decay amplitudes of the one-gluon exchange and the con- 
finement interactions should be added coherently. Therefore, the current- 



current interactions can be written in the generic form as [15 



H, 



1 



d\d\ r{x)K{\x - y\)r{ 



(9) 



The current J" in Eq. ([9]) is assumed to be a color octet. The currents, J, 
with the color dependence A"/2 factored out and the kernels, K{r), for the 
interactions are 



Currents 



ip{x)I^jj{x) 
%l>{x) 7° il){x) 



Kernels 



K{r) 



Scalar Lorentz current. 
Time-like vector Lorentz current. 
Space-like vector Lorentz current. 

(10) 



Scalar Confinement, 



-4a, [-a,(l - + A] 

+4(1 — a,) [— ac(l — e"^"^) + A] Static Vector Confinement, 
-4(1 - as) [- 



I eta 

r 



'1 — e ^'''^) + A] Transverse Vector Confinement, 
Coulomb OGE, 
Transverse OGE. 

(11) 
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Following Ref . [15| , we refer to this general type of interaction as a JKJ 
and to the specific cases considered here as sKs, j^Kj^ and j'^Kj^ interac- 
tions. Details of the resulting matrix elements for different cases are given 
in 



Appendix A 



The diagrams that contribute to the strong decay A B + C are shown 
in Fig. [TJ There are two coming from the quark line, dig and d2q. The 
difference between them is the rearrangement of the quarks and antiquarks 
in the final mesons. The other two diagrams are referred to the antiquark 
line, dig and ^25- 

The total width is the sum over the partial widths characterized by the 
total spin, Jbc, and the relative angular momentum, I, of the final mesons 
B and C 

^A^BC= ^A^BcjJBcJ), (12) 

Jbc^I 

where 

TA^BciJBcl) =2n J dko5{EA - EBc)\MA^Bc{k^)? (13) 



and M.A^Bc{ko) is calculated following Appendix A 
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Using relativistic phase-space we arrive at 



^a-^bc{Jbc, 



2n^^\MA^Bciko)\' 
niAko 



where 



kn 



2mA 



(14) 



(15) 



is the on-shell relative momentum of mesons B and C. 

We will compare our results with the widely use ^Pq decay model which is 
a particular case of Eq. f lTTj) with only a constant scalar term 5*0 = ^^mq^y. 
A complete development of this model can be found in Ref. 



4. Results 

From an experimental point of view there are a few data in the open- 
charm decays of the 1 cc resonances. The main experimental data are the 
resonance parameters, mass and total decay width, of the excited ip states. 

Tables [3] and m show the strong decay widths predicted by the microscopic 
model for the 1 cc states established in Table [2] compared with the exper- 
imental data and the ^Pq results. The parameter 7 of the ^Pq model has 



been fitted in Ref. [29|. The notation D1D2 includes the D1D2 and -D1-D2 
combination of well defined CP quantum numbers. For the kinematics we 
use experimental masses whenever they are available. 

One can see that the total decay widths predicted by the microscopic 
model are in general lower than the experimental ones, whereas those pre- 
dicted by the ^Pq model reproduce the data in a better way. However, it is 
worth to notice that the correct order of magnitude of the strong decays is 
given by the microscopic model with no free parameter in contrast with the 
one parameter of the ^Pq model. 

When we go up through the spectrum, the states are more and more wide 
and the total widths for S and D-waves are larger in both decay models, 
always D-wave widths are smaller. 

In order to disentangle the contribution of the different quark-quark po- 
tential pieces, we compare in Table [5] the results of the full model with the 
ones taken into account only the time-like component of the confinement 



potential. These last results can be compared with those of Ref. [lj| which 



includes the same pieces of the current although in a slightly different model. 
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It seems that results including only time-like vector confinement compo- 
nent are better than those of the full model. However, if one looks to the 
ratios between the different decay channels (Table |6]) none of the models are 
able to reproduce the experimental data. This fact suggest that the dynamics 
of the charmonium strong decays is far from be a simple process and more 
Fock components of the wave function can be involved in the decay [30] . 

5. Conclusions 

Microscopic models of meson strong decays into two mesons depend on 
the transition Hamiltonian which drives the decay mechanism. We have de- 
veloped a model in which the full Hamiltonian that determines the spectrum 
is used for the decay. In general, the obtained total decay widths are lower 
than the experimental data although the order of magnitude is reproduced 
without any free parameter. 

It seems that, considering only the confinement time-like components, 
the agreement with the experimental data is improved. This fact seems to 
be in line with the conclusions of Ref. |[21j]. The authors stated that the 
Dirac structure of confinement should be of a time-like nature which dy- 
namically generates an effective scalar interaction as required by the hadron 
spectroscopy. However, fine details of the charmonium decays, like the ratios 
between different decay channels, are not reproduced by any model which 
suggests that the strong decays into charmed mesons is still an open prob- 
lem. 
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Appendix A. Matrix elements in the microscopic model 

Appendix A.l. Transition operator 

If one considers only the contributions in which a quark- ant iquark pair 
is created, the interaction Hamiltonian in Eq. ([9]) reduces to the following 
transition operator 

dW-K{\x-y\) f d^r^p^^^ J^^^Lyyyy 



+ al{Pi)arAP2)al,{pMM [urAPi)Tu,AP2)] [urM)^^rAP^)] (T''^^-^^-' (T^^^^^^-y 

+ al{Pl)hl{p2)hrMK{P^) [UrAPl)TVrM)] [VrM)^^rAP^)] e-'^^^^^'^ C^'^^-^^'y 

+al{pi)hl{P2)al{p,)a,M) [urAPi)TvrAm [urM)^^,^)] e-(^^^+^^^)-V^(^^3-P4)-.^] 

(A.l) 

where the first term is equal to the third one. This can be seen exchanging 
the X and y variables in the first term and then, changing 1^3 and 2^4 
particles taking into account the anti-commutation rules of the creation and 
destruction operators to arrive to the third term. This is possible because the 
kernel depends on x and ^ as \x — y\. The same occurs with the second and 
fourth terms. Therefore we have a factor two and we can write the transition 
operator as 

rr /"j3 r3 r^/i- -ia f d%d%d%d% I mim2mzmi \- 
1=1 dxdy K[\x — y\) / — — w— — - — - — -— > 

J ^ ^1 (27r)6 \j E,^,Ep,Ep,Ep, ^^^^^^^^ 

+alM)hl{p2)hrApMM) [UrAPl)TVrAP2)] [t^.a (P^rt^.^ (P4)] e-^^^^^^^^^)' V^^^^^'"^^^)'^^] , 

(A.2) 

where the first and second terms refer to the qq pair creation from the quark 
line and from the antiquark line, respectively. The diagram representation of 
these two terms can be seen in Fig. [H diagrams diq and diq. For illustration 
we build the result from the diagram dig, the transition operator is 



j3 j3 T^f\^ -i\ f d^id%d%d% riuYm^m^m^ \- 
dxdyK(\x — y\) / -, — — J > 

[ al^{p,)bl{p2)al{ps)arM) [urAPi)^VrM)] [nrM)^^^rM)] e'^^^^^^^^^' 6-^^^^-^^^^^ 

(A.3) 
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The calculation of the diagram dig can be followed from that of the diagram 
dig. If the initial meson is formed by a quark and an antiquark with equal 
masses, the contribution of both diagrams to the decay rate is the same and 
they contribute constructively. 

Now, we can integrate in x and y 

T^J d%d%d%d%K{\Q\)5^'\pi+p2+P3-p4] 



mim2m3m4 

-E'pi -E'p2 -^P3 ■^P4 



X] [41(^1)^^2(^2)43(^3)0^4(^4) Ki(pi)rvr2(P2)] K3(P3)rMr4(P4)]] , 

ri,r2,r3,r4 

(A.4) 

where Q = pi + P2 = p^ — Ps is the momentum transferred, f^dQI) is the 
Fourier transform of the kernel K{r) and the 5-function implies the momen- 
tum conservation. 

Appendix A. 2. Transition amplitude 

We are interested on the transition amplitude for the reaction {a 13) a — )■ 
{^^)b + {^p)c- In the center-of-mass reference system of meson A one has 
Ka = Kq = and the matrix element factorizes as follow 

{BC\T\A) = S^^\Ko)Ma^bc. (A.5) 
The initial state in second quantization is 

1^) = J d'p^d%S^'\KA - pA)MPA)ai{Pa)bl{p0) |0) , (A.6) 

where a {(3) are the spin, flavor and color quantum numbers of the quark 
(antiquark). The wave function 0a (pa) denotes a meson A in a color singlet 
with an isospin I a with projection M/^, a total angular momentum J a with 
projection Ma, J a is the coupling of angular momentum La and spin Sa- 
The Pa and Pj3 are the momentum of quark and antiquark, respectively. The 
Pa and pa are the total and relative momentum of the {af3) quark-antiquark 
pair within the meson A. The final state is more complicated than the initial 
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one because it is a two-meson state. It can be written as 



\BC) = \ = / d?KB(fKc V {JBcMBclm\JTMT)6^^\K~K^)5{k-h) 

y^ + ^BcJ ^^^^^ 



^""^^^ {JbMbJcMc\JbcMbc) {IbMiJcMiJIaMi^) 



k 

MbMcMibMic 



d^psd^pJ^Pxd%S^^\KB - Pb)S^^\Kc - Pc) 
(pB{pB)(j)c{pc)al{p5)bl{pe)a{{px)bl{pf,) |0) , 



(A.7) 



where we have followed the notation of meson A for the mesons B and C. 
We assume that the final state of mesons B and C is a spherical wave with 
angular momentum /. The relative and total momentum of mesons B and 
C are ko and Kq. The total spin Jbc is obtained coupling the total angular 
momentum of mesons B and C, and Jt is the coupling of Jbc and /. 

The diagrams that contribute to the reaction and are allowed by the 
transition operator are shown in Fig. [H Two of them are coming from the 
quark line, dig and d2q, and take into account the different rearrangement of 
the quarks and antiquarks in the final mesons. The other two diagrams are 
referred to the antiquark line, dig and d2g. We have different cases: 

• Case in which a = fi = (3. The two diagrams, dig and d2g, contribute to 
the decay amplitude. The contribution of diagram dig is Mji^bc and 
the contribution from diagram d2q can be calculated from the amplitude 
of the dig diagram changing meson B and C (Ma^cb), so the total 
amplitude is given by 

Ma-.bc = Ma^bc + {-IY^+'^-'^+'^+'^'-'^^+'Ma^cb. (A.8) 

• Other case. Only one of the two diagrams contribute to the amplitude 

M A-^BC — Ma^bc- (A. 9) 

If the quark and antiquark in the original meson are the same, then the 
contribution of diagram dig {d2g) is equal to the diagram dig {d2g) and both 
contribute constructively. In other case they have to be calculated separately. 
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When the initial A meson has definite C-parity we have to use final states 
with definite C-parity. If CB = C the state has definite C-parity and the 
amplitude is given by the above rules. If CB ^ C then the appropriate 
C-parity combination has to be taken and this gives a factor \/2 in the 
amplitude (or the amplitude cancels for the wrong C-parity). 

For illustration we build the result from the diagram diq (Ma^bc)- The 
amplitude is a product of a Fermi signature phase, a color factor, a flavor 
factor and a spin-space overlap integral 

■^A^BC -^signature ^ -^;olor ^ -5ja,vor ^ -^spin— space ■ (A. 10) 

Appendix A. 2.1. Fermi signature phase 

The Fermi signature can be read off from the diagram as the number of 
line crossings because it arises from the ordering of the quark and antiquark 
operators. In the case of dig diagram we have 

/signature = (-1)' = -1- (A.ll) 

Appendix A. 2. 2. Color factor 

As all mesons are color singlet qq states the color term is given by 



r - —\^Tr 
^color 3 / ^ -"-^ 

32 „ 



3i- 



(A.12) 



Appendix A. 2. 3. Flavor factor 
For the flavor sector we have 



^flavor - (-1)*«+*''+^V(2/b + 1)(2/C + 1) {J^ \l J^}, (A.13) 

where is the isospin of a given quark or antiquark ^. Note that the isospin 
operator in the creation vertex is uu-\- dd-\- ss. 
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Appendix A. 2. 4- Spin- space factor 

The spin-space overlap integral for the diagram diq. 1 •<->■ 2 •<->■ i/, 3 ■<->■ 5' 
and 4 -H- q;', reads as follow 

-^spin— space ^ ^= 



d'^KBd'^Kcd%d%d'^Pxd\d'^Pi,d%d%,d'^p^,d'^Pad% f^^^^^^ 

k 

<pBiPB)(i)c{Pc)(i)A{PA) K{\pi^+p^\)d^^\pi^+p^+ps' -Pa') 

J2 ^c'a^^'^Pc^' -Pa)Ss5'S^'\p5-P5')Se.S^'\pe-p.)Sx,5^'\px-p^) 

a' ,5' ,(1,1/ 

^pp5^^\Pp-Pf}) [u^i{.P^,)Tv„{p^)\ [u5>{p5>)TUa'{pa')] ■ 

(A. 14) 

Now using some 5-functions in momentum and spin of quarks (antiquarks), 
we can simplify the above expression 

-^spin— spcLce ,-. = 

{JBcMBclm\JTMT){JBMBJcMc\JBcMBc) 
[ d'KBd'Kcd%d%d'p,d'p^d'p^d%.l^^^^^ 

J V ^Pp^^Pu^PS^Pa 

k 

(I>b(pb)(I)c(pc)(I>a(pa) K(\p^+p^\)5^^\ps - (pa -Pi^- Pi.)) 

5pp5^^\Pp-Pp) [u^i.{Pn)^V„{p,,)\ [u5{p5)Vu^{po)\ . 

(A.15) 
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The nonrelativistic reduction of Eq. (lA.lSp without specifying the JKJ decay 
model is 

1 



spin— space n. j 

V 1 + Obc 



{JBcMBclm\JTMT) {JbMbJcMcUbcMbc) 

m,MBC,MB,Mc 



k 



(fKBd^Kcd%(fppd^Pf,d%d^Pad^Pp6'-^\K - Ko)6{k - ko 

S^^\Kb - Pb)S^'\Kc - Pc)s'-'\Pa)Mpb)Mpc)Mpa) 

K{\P^ + Pu\)S^^\p5 - (Pa - - Pu))SppS^^\pp - Pp) 

hm [up{pf,)Tvy{py)] hm [us{p&)Tua{pa)] ■ 

(A.16) 

We require spin matrix elements which involve the nonrelativistic 0{p/m) 
matrix elements of Dirac bilinear s with F = J, 7°, 7 and Pauli spin matrix 
elements. These are 

{.Pu-p^,) ■ 



lim [u^,{pp) Iv^{p^,)] 

t)/c— i-O 



2m 



2m,, 



-V25 



TT 



3^1 



Pli-Pu 



22 



lim [up{pp)-i^v^{p^)\ 



2m, 



(A.17) 



2m„, 



TC 



yi{Pu+Pt,) 



1 1 
22 



lim [up{pp,)^v^{p^)] = 



and 



lim [Uf,{pf,) I u^{p,y)] =6/,^, 
lim [u^,{p^)l'^u„{p^)] = 

v/c— i-O 



^limJMP.)7-.(p.)] = ^ 



where we have used the relation 

fi-((T) = -V3[ri®(a)]o = y6 



1 1 
22 



(A.18) 
(A.19) 
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Then, the expression for the different contributions are 
• sKs interaction 

-1 1 



1!. 



spm-space ^^y+S^ 2m, 



V2^ / d^KBd^Kcd^Pad^ppSpi^d^pj.d^ps 



S'-^^ {PS - {Pa -Pfi- Pv)) K{\p^ +p,\) 

{ [[(l>B{PB){SaSi.)SB]JB[(l>c{Pc){Sij.Sp)Sc]Jc]JBCyi{k) Jt\ 



^A{pA){SaSi3)SA] J A 



yi 



Pix-Pu 



Ja). 
(A.20) 



j^Kj*^ interactions 



spm-space ^/^qi^ 2m, 



\/2% / d^KBd^Kcd^pad%d^p^d^p^d% 



6^'\K - Ko)S^'HKb - Pb)S^'\Kc - Pc)S^'\Pa)^-^^j^ 

5(3) {P5 - {Pa -Pf,- Pu)) K{\p^ + p,\) 

{ [[(f)B{pB){SaSu)SB] Jb [0c(pb) (Sm*/3)'S'c] M JBcYl{k) 

I [[0A(PA)(saS/3)-S'A] Ja [yi {Pix+Pv) (s^s^)l] 0] Ja) ■ 



Jt\ 
(A.21) 



j^Kj'^ interactions 



1 



spm-space ^/Y+I^ 



J2 {JBcMBclm\JTMT) {JbMbJcMc\JbcMbc) 

m,MBc,MB,Mc 

YimCk) 



J d^KBd^Kcd^psd^ppd^Pf,d^p,d%d^p^ 5^^\K - Ko)S{k - ko 

S^^\Kb - Pb)5^'\Kc Pc)^^''\Pa)Mpb)Mpc)Mpa) 
K{% + p^)5^'^\p5 - {Pa -p; - Pu))5ppS^^\pp - pp) 



k 



lim [u^{pi^)Yv,{p,)] 5ij - 



) lim [us{psh^Uc,{pa)] . 

(A.22) 
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The procedure followed to solve the above spin-space overlap integrals is 
similar to that of Ref. 28 . 
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Meson 


State 


Channel 








^Mic. 


^(3770) 


l^Di 


D+D- 


11.34 


42.8 


8.03 


42.3 








15.13 


57.2 


10.94 


57.7 






DD 


26.47 


100 


18.97 


100 


27 c _|_ f) Q 




total 


26 47 




18 97 




^(4040) 


3=^5i 


DD 


4.61 


4.1 


10.17 


26.0 






DD* 


22.23 


20.0 


18.75 


47.9 






D*D* 


82.35 


74.0 


9.06 


23.2 






D,D, 


2.08 


1.9 


1.14 


2.9 


80 ± 10 




total 


111.27 




39.12 




^^(4160) 


23^1 


DD 


22.82 


19.7 


17.03 


52.1 






DD* 


2.22 


1.9 


7.38 


22.6 






D*D* 


83.73 


72.2 


5.28 


16.2 






D,Ds 


0.24 


0.2 


2.61 


7.9 






DsD* 


6.94 


6.0 


0.40 


1.2 


103 ±8 




total 


115.95 




32.70 






^ u\ 


n n 




7 n 


0. / o 


^ P. 

0.0 






DD* 


8.19 


7.2 


29.81 


29.2 






D*D* 


8.87 


7.8 


46.46 


45.5 






DD^ 


54.51 


47.8 


2.18 


2.1 






DD[ 


4.29 


3.8 


12.02 


11.7 






DDI 


27.17 


23.8 


0.56 


0.6 






D,D, 


0.07 


0.1 


1.86 


1.8 






DsDl 


1.90 


1.7 


3.36 


3.3 






d:di 


0.91 


0.8 


0.17 


0.2 


74 ± 15 ± 10 




total 


113.92 




102.15 




V'(4415) 


33^1 


DD 


15.11 


9.5 


7.93 


18.5 






DD* 


5.82 


3.7 


6.66 


15.6 






D*D* 


32.56 


20.5 


7.23 


16.9 






DDi 


64.77 


40.7 


6.06 


14.2 






DD[ 


6.92 


4.4 


2.12 


5.0 






DDI 


23.60 


14.8 


1.82 


4.3 






D*D* 


7.12 


4.5 


2.39 


5.6 






DsDs 


0.31 


0.2 


2.22 


5.2 






DsD* 


0.68 


0.4 


1.09 


2.5 






DlDl 


2.13 


1.3 


5.20 


12.2 


62 ±20 




total 


159.01 




42.72 





Table 3: Open-flavor strong decay widths, in MeV, and branchings, in %, of -0 states. 
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Meson State 
X{A630) 5^ 



Q9+40+10 

X(4660) 43Di 



48 ± 15 ± 3 



Channel 


Tap,, 








DD 


6.62 


3.2 


1.44 


0.8 


DD* 


26.23 


12.7 


15.82 


8.4 


DW 


15.57 


7.5 


30.40 


16.2 


DDi 


2.88 


1.4 


18.70 


9.9 


DD[ 


4.52 


2.2 


2.58 


1.4 


DDi 


0.00 


0.0 


21.14 


11.2 


D*D* 


6.97 


3.4 


10.10 


5.4 


D*Di 


39.21 


19.0 


22.47 


11.9 


D*D[ 


14.35 


7.0 


26.24 


13.9 


D*Dl 


80.47 


39.0 


18.28 


9.7 


DsD, 


0.92 


0.4 


1.28 


0.7 


DsD* 


0.30 


0.1 


6.70 


3.6 


D*D* 


1.14 


0.6 


6.34 


3.4 


^ S-'-^SL 


2.82 


1.4 


0.92 


0.5 


D,D\ 


0.79 


0.4 


0.03 


0.0 


D.D% 


0.19 


0.1 


0.22 


0.1 




2.76 


1.3 


1.30 


0.7 


si 


0.14 


0.1 


3.74 


2.0 


^s^sl 


0.26 


0.1 


0.29 


0.1 


-^sO sO 


0.22 


0.1 


0.23 


0.1 


total 


206.37 




188.22 




DD 


10.92 


8.1 


3.21 


2.3 


DD* 


7.55 


5.6 


4.10 


2.9 


D*D* 


38.04 


28.2 


2.67 


1.9 


DDi 


2.41 


1.8 


20.51 


14.4 


DD[ 


0.51 


0.4 


2.62 


1.8 


DDI 


0.00 


0.0 


6.75 


4.8 


D*Dl 


3.44 


2.5 


0.71 


0.5 


D*Di 


34.83 


25.8 


10.89 


7.7 


D*D[ 


6.98 


5.1 


2.96 


2.1 


D*D* 


21.92 


16.2 


77.52 


54.5 


DsD, 


0.96 


0.7 


1.46 


1.0 


DsDl 


0.00 


0.0 


1.35 


0.9 


DtD* 


0.33 


0.2 


4.28 


3.0 


DsDsi 


3.63 


2.7 


0.0 


0.0 




1.09 


0.8 


0.62 


0.4 


DsDt, 


0.08 


0.1 


0.07 


0.1 


did:. 


1.18 


0.9 


0.43 


0.3 


DID,, 


0.48 


0.4 


0.93 


0.6 


DID',, 


0.17 


0.1 


0.37 


0.3 


r>* D* 


0.53 


0.4 


0.74 


0.5 


total 


135.06 




142.19 





Table 4: Open-flavor strong decay widths, in MeV, and branchings, in %, of tp states 
(Continuation). 
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Decay 


Ref. 


14 


fKf 


Mic. 


Exp. [25] 


^(3770) ^ DD 


20. 


1 


29.8 


19.0 


25.6 ±3.4 


total 


20. 


1 


29.8 


19.0 


27.5 ±0.9 


^^(4040) ^ DD 


0.1 


1.4 


10.2 




^'(4040) DD* 


33. 





25.2 


18.7 




^'(4040) D*D* 


33. 





35.0 


9.1 




^'(4040) DsDs 


8.0 


0.3 


1.1 




total 


74. 





61.9 


39.1 


80± 10 


^{4160) DD 


3.2 


25.0 


17.0 




V'(4160) DD* 


6.9 


0.5 


7.4 




^p{AlQO) D*D* 


41. 


9 


21.3 


5.3 




^'(4160) -> DsDs 


5.6 


0.03 


2.6 




i){A160) -> D,D* 


11. 





0.6 


0.4 




total 


69. 


2 


47.4 


32.7 


103 ±8 



Table 5: Open-flavor strong decay widths, in MeV, of ip states reported in Ref. fl3] and 
our decay rates taking into account the static vector contribution or the full model. 



State 


Ratio 


fKf 


Mic. 


3p 


Ref. [14] 


Measured [25] 


^'(4040) 


DD/DD* 


0.06 


0.54 


0.21 


0.003 


0.24 ±0.05 ±0.12 




D*D*/DD* 


1.39 


0.48 


3.70 


1.0 


0.18 ±0.14 ±0.03 


V'(4160) 


DD/D*D* 


1.17 


3.23 


0.27 


0.076 


0.02 ±0.03 ±0.02 




DD*/D*D* 


0.02 


1.40 


0.03 


0.16 


0.34 ±0.14 ±0.05 


ip{U15) 


DD/D*D* 


1.54 


1.10 


0.46 




0.14 ±0.12 ±0.03 




DD*/D*D* 


0.28 


0.92 


0.18 




0.17 ±0.25 ±0.03 



Table 6: Open- flavor strong ratios of ip states predicted by different decay models and 
their comparison with the experimental data. 
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